Thermal valence- bond- solid transition of quantum spins in two dimensions 
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We study the S = 1/2 Heisenberg (J) model on the two-dimensional square lattice in the presence 
of additional higher-order spin interactions (Q) which lead to a valence-bond-solid (VBS) ground 
state. Using quantum Monte Carlo simulations, we analyze the thermal VBS transition. We find 
continuously varying exponents, with the correlation-length exponent v close to the Ising value for 
large Q/J and diverging when Q/J approaches the quantum-critical point (the critical temperature 
T c — » 0). This is in accord with the theory of deconfined quantum-critical points, which predicts 
that the transition should approach a Kosterlitz-Thouless (KT) fixed point when T c — >• + (while 
the transition versus Q/ J for T = is in a different class) . We find explicit evidence for KT physics 
by studying the emergence of U(l) symmetry of the order parameter at T = T c when T c — » 0. 

PACS numbers: 75.10.Kt, 75.10.Jm, 75.40.Mg, 75.40.Cx 
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The S = 1/2 Heisenberg model on the two-dimensional 
(2D) square lattice can host a quantum phase transition 
between the standard Neel antiferromagnet (AFM) and 
a valence-bond-solid (VBS) ground state when other in- 
teractions are added [1]. This transition between two 
different ordered quantum states has been the subject 
of a large body of work for more than 20 years Q. In 
the J-Q model [3|, the pair exchange J is supplemented 
by a product of two or more singlet-projectors on adja- 
cent links of the lattice, with strength Q. For sufficiently 
large Q/J, the correlated singlets destroy the Neel order 
existing for small Q/J, leading to the VBS crystalliza- 
tion of ordered singlets. Unlike geometrically frustrated 
systems, on which searches for VBS states and the AFM- 
VBS transition were focused for a long time 0-Q ? the J~ 
Q model is amenable to large-scale quantum Monte Carlo 
(QMC) simulations ^ and its AFM- VBS transition has 
been studied extensively [H, I9l-[l8|. Many results indi- 
cate that the model realizes the unusual ( "non-Landau" ) 
deconfined quantum-critical (DQC) point proposed by 
Senthil et al. [HI, [lo|, where the two order parameters 
both arise out of of emergent spin- 1/2 degrees of freedom 
(spinons), which at criticality are described by a gauge- 
field theory; the non-compact CP 1 model. Other, less 
exotic scenarios within the standard Landau-Ginzburg- 
Wilson framework for phase transitions have also been 
put forward [111, E0, E2 , however. 

The putative DQC points are manifestations of inter- 
esting quantum effects, due to Berry phases and emergent 
topological conservation laws [lo|, HH, that potentially 
are at play in many strongly-correlated quantum-matter 
systems. Being amenable to large-scale unbiased QMC 
simulations, further studies of the J-Q class of models 
offer opportunities to examine the DQC proposal in de- 
tail from various angles. Here we present results for the 
ordering transition of the VBS at finite temperature, dis- 
cussing its universality, relationship to conformal field 
theory (CFT), and insights gained into the emergent U(l) 
symmetry [20| associated with the DQC point when ap- 



proached at finite temperature. 

Universality of the VBS transition — The square-lattice 
columnar VBS obtaining with the standard J-Q model 
breaks Z4 symmetry and, thus, it should also exist at 
finite temperature (T > 0). Thermal 2D ^-breaking 
transitions normally do not have fixed critical exponents, 
but belong to a universality class of CFTs with charge 
c = 1 exhibiting continuously varying exponents (as a 
function of model parameters) [24|, [25| . Realizations of 
these transitions include the standard XY model with 
a field fecos(4^) for all spins i (angles Qj) 26|, 27], the 



Ashkin- Teller model [28|, |29| , and the Ising model with 



nearest- and next-nearest neighbor interactions (the J\- 



J2 model) [30|,[3l|. The deformed XY model has a critical 



line connecting Ising and Kosterlitz-Thouless (KT) fixed 
points [H, [33|, while the critical lines of the AT and J\- 
J2 models connect Ising and 4-state Potts points. It is 
then intersting to ask if any of these scenarios are re- 
alized in the T > paramegnet-VBS transitions of the 
J-Q model. In this Letter we present strong evidence 
for universality corresponding to the Ising-KT critical 
line, with the KT transition obtaining in the limit when 
Q/J approaches its quantum-critical value and the crit- 
ical temperature T c 0. This is in agreement with the 
DQC theory and its U(l) gauge-field description, where 
the nature of the VBS state is dictated by a dangerously 
irrelevant operator [1, [H, [lo| , which implies that the VBS 
fluctuations should cross over from Z4 to U(l) symmetric 
as the DQC point is approached, which in fact has been 
observed in ground state studies of the VBS fluctuations 
of J-Q models d, Qj], 12 1- We here show explicitly that 
this also applies to the T > critcal line when T c — >> 0. 

The T > VBS transition was previously studied by 
Tsukamoto, Harada and Kawashima [34j |. who carried 
out QMC simulations of the J-Q2 version of the J-Q 
model, where the Q2 interaction is one of products of 
two singlet projectors. The results were puzzling, with 
significant deviations from the "weak universality" sce- 
nario applying to the transitions discussed above, where 
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the critical correlation- function exponent rj = 1/4 is con- 
stant (while other exponents depend on system details). 
Instead, 77 « 0.5 was obtained [34|. Here we consider 
the J-Q3 model jl2j |. where the Q3 term consists of three 
bond-singlet projectors (forming columns on three adja- 
cent lattuce links). This model has a much more robust 
T = VBS for large Q3, while the VBS state of the 
J-Q2 model is near-critical even for Q2/J —> 00. With 
the J-Qs model we can systematically study the T > 
transition both far away from the DQC point and close 
to it. We find consistency with 77 = 1/4 to high precision, 
and also point out that cross-over behavior related to the 
DQC criticality exactly at T = makes it difficult to re- 
liably extract the exponents when T c is low. We believe 
that this behavior affected the previous study of 77. 

Model and methods — We next discuss the QMC calcu- 
lations and data analysis on which we base our conclu- 
sions. The J-Qs Hamiltonian is defined as 
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where is a nearest-neighbor bond-singlet projector; 
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here on the square lattice with L 2 sites. We define the 
coupling ratio q = Qs/J- The quantum-critical point 
separating the AFM and VBS states is q c = 1.500(2) (T3|. 
We here use the stochastic series expansion (SSE) QMC 
method with loop updates [35l437j to compute several 
quantities useful for extracting the critical temperature 
and exponents of the VBS transition for q > q c . 

There are various ways to define the VBS correlation 
length. For computational convenience we here use a 
definition based on the J-term (bond) susceptibility, 



Xbi ,6 2 



dT(P b2 (T)P bl (0)), 



(3) 



where P5 is a singlet projector as in (j2j), with b denoting a 
bond connecting sites ib,jb- These susceptibilities can be 
computed easily with the SSE method, because the bond 
operators are terms of the Hamiltonian and, thus, appear 
in the sampled operator sequences. With n(b) denoting 
the number of J-operators on bond b in the sequence, the 
susceptibility is given by [HI 



Xb u b 2 = {n{bi)n(b 2 ) - S blib2 n(bi)) / f3 . 



(4) 



This estimator works well as long as q is not too large. 
When q > 10 the measurements become noisy due to the 
low density of bond operators, but for our purposes here 
this is not a problem. 

To detect columnar VBS order, we consider the bonds 
b\ and 62 oriented in the same (x or y) lattice direction 
and denote by x a ( Y )i a = x iVi the spatially averaged 
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FIG. 1. (Color online) Extraction of T c for system at q = 5. 
Shown in (a) are, in order of higher to lower curves on the left 
side, results for £i/L versus T for system sizes L = 96, 48, 24, 
and 12. Crossing points giving T C (L) estimates are shown 
in (b), using both £1 and £2 with sizee pairs (L,2L). The 
data were fit to the form T C (L) = T c (oo) + a/L w in the range 
1/L e [0,0.08] (£1) and [0,0.06] (&), yelding T c = 0.249(3) in 
the case of xi- For the £2 fit, T c (oo) = 0.249 was fixed. 

distance-dependent susceptibility. The VBS susceptibil- 
ity Xvbs i s ^ ne q = (^j 0) Fourier transform of x x ( r ) (and 
analogously for y). Because the columnar VBS breaks 
the lattice rotational symmetry, we can define two cor- 
relation lengths. Using the x susceptibility and defining 
q = (tt,0), qi = (7r + 27r/L,0) and q 2 = (7r,27r/L) we 
have the correlation lengths parallel and perpendicular 
to the x-oriented bonds for an L x L lattice; 



tx = L_ / XvBs(qo) 

2^Yxv-Bs(qi) 



tx-lL / XvBs(qo) 1 

' 42 2^VxvBs(q2)" 



(5) 

and analogously for y. Average valuess of x, y quantities 
are denoted in the following without superscript. 

Critical temperature — To illustrate how the critical 
VBS temperature T c is determined, Fig.QJa) shows 
versus T at q = 5 for several system sizes. According 
to standard finite-size scaling theory [39|, £i/L for differ- 
ent L should cross at T c when L 00. Due to expected 
scaling corrections, the crossing point T c (Li, L2) between 
two system sizes, which we here take as L and 2L, drifts 
slowly with L and converges as the system size increases. 
We use the crossing point for both £1 and £2 to extract 
T c and check the consistency of the two results. 

Fig.QJb) shows two sets of T C (L) point obtaied from £1 
and £2- Both curves can be fitted with the form T C (L) = 
T c (oo) + a/L w but the parameters are different. The 
two curves appoach T c from different directions. The £1 
data have large deviations from the fitted function only 
for small systems (L < 12), while £2 shows corrections 
extending up to larger systems and the size dependence 
is non-monotonic. In spite of the different behaviors, 
the data extrapolate consistently to a common T c in the 
thermodynamic limit. To demonstrate this, we show in 
Fig. Hp)) a fit to the £1 data, which gives T c = 0.249(3). 
(which has a smaller statistical error than the value from 
£2). We also show a fit to the £2 data, where the T c (oo) 
value is fixed at the result based on £1. 

T c values for several other q points were extracted in 
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FIG. 2. (Color online) (a) The critical temperature extracted 
from £i/T (open circles). Also shown are results (solid circles) 
where the VBS susceptibility exhibits the best scaling behav- 
ior when 7 = 7/4 is fixed, (b) The exponent v versus q. The 
vertical dashed lines in both panels mark the quantum-critical 
ratio q c [13] • The curves are guides to the eye. 

the same way, making sure that £i and £2 data extrap- 
olate consistently but using only the £1 results (which 
always have smaller errors) for further analysis. This pro- 
cedure becomes increasingly challenging as the quantum- 
critical point q c is approached and T c 0. The correc- 
tions to the asymptotic form became more profound and 
larger systems have to be used. In addition, the SSE cal- 
culations become more time-consuming, since L ^> 1/T 
is required for the simulated effective classical system to 
be firmly in the 2D limit. The largest system simulated 
was L = 192atg = 5/3. Results for T c are shown versus 
the coupling ratio in Fig. [2ja). 

Critical exponents — we next present an analysis of the 
scaling behavior of the VBS susceptibility, which exactly 
at T c should follow the form 



XvBS^c)-^, 



(6) 



where 7/^ = 2 — rj. Here we can use the value of T c 
extracted above from the correlation length scaling. Al- 
ternatively, we can adjust the temperature until the best 
power-law scaling is obtained. If sufficiently large sys- 
tem sizes are used the two methods should of course de- 
liver consistent results. This is indeed the case, as shown 
in Fig. Efa). An example of the best power-law scal- 
ing is shown for the system with q = 5 in Fig. [3^a). 
Here the corrections to scaling appear to be very small 
(i.e., a straight line can be well fitted on the log-log scale 
even when systems as small as L = 10 are included) and 
the temperature, T = 0.253, is only about one error bar 
off the T c value extracted from £i/L. A series of fits 
with a bootstrap analysis to estimate the errors yielded 
j/v = 1.750(1), corresponding n = 0.250(1). Thus, we 
find complete consistency, to rather high precision, with 
the most natural expectation of n = 1/4. We obtain 
similar results for all values of Qs/J studied. 

Fig. [3fb) demonstrates a different way to analyze the 
susceptibility and test the assumption n = 1/4, by graph- 
ing Xvbs£ -7 ^ 4 versus T is for different system sizes. All 
curves cross essentially at the same point, which confirms 
the scaling power 7/z/ = 7/4 in Eq. (|6]). The remarkable 
absence of drift in the crossing points of Xvbs^ -7 ^ 4 (in 
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FIG. 3. (Color online) (a) Scaling behavior of the critical 
VBS susceptibility for systems at q = 5. Here T was ad- 
justed to give the best linear scaling on the log- log plot, giv- 
ing 7/f = 1.750(1). (b) The size-scaled susceptibility under 
the assumption n = 1/4 versus T for several system sizes. 
The crossing point is consistent with T c extracted from the 
correlation length. 
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FIG. 4. (Color online) Data collapse of the VBS susceptibility 
for system s at q — 10/3. The inset shows data for L — 
80, 96, 112 in the range tL 1/u e [-0.5, 3] for which the fitting 
procedure was carried out. The main part shows data in a 
larger window and including also smaller systems. The fit 
yelded v = 1.70(5). 

contrast to the significant drift found for the normalized 
correlations lengths) makes this quantity a perfect candi- 
date for carrying out a finite-size data collapse to extract 
correlation length exponent ^, which we consider next. 

Shown in Fig. |4] are data sets for system sizes L = 48 
to 112 at q = 10/3, graphed versus tL 1 /" , where t is the 
reduced temperature, t = (T — T c )/T c , and the critical 
temperature was determined in the manner above to be 
T c = 0.217. The correlation lengt v was adjusted to 
give the best data collapse, as measured with respect 
to a polynomial fitted simultaneously to all data points 
for L = 80,96,112 in the range tL 1 ^ e [-0.5,3]. A 
zoom-in on this window is shown in the inset. The fit 
was restricted to the larger sizes in order to minimize the 
effects of neglected scaling corrections, and the window of 
tL x l v values was chosen to obtain a statistically sound fit. 
This procedure along with an analysis of the statistical 
errors gave v = 1.70(5). When q is tuned towards g c , 
larger system sizes are required to achieve good collapse 
due to more pronounced scaling corrections, as already 
mentioned above. As an example, at q = 5/3, we used 
system sizes L = 112, 128, 160, 192. 

All our results for T c and v versus q are shown in Fig. [2] 



T c clearly decreases when q approaches q c and v grows 
rapidly, changing from 1.065(5) at q = 10 to 2.7(1) at 
q = 5/3. The behavior suggests that v diverges when 
q — >• <2c, which would mean that the critical line corre- 
sponds to the c = 1 Ising-KT scenario, with the KT 
universality applying in the limit q — >• q^ and 2D Ising 
universality (y = 1) applying in the extreme limit far 
from the quantum- critical point (which cannot strictly 
be achieved within the J-Qs model, but v is already close 
to the Ising value for q = 10; the largest q studied here). 
This scenario is also supported by the fact that there is 
no specific-heat peak at T c , i.e., the exponent a < 0. 

Emergent U(l) symmetry — The changing critical ex- 
ponents are related to an evolution of the critical VBS 
fluctuations. We investigate these by following the dis- 
tribution of the components (D x ,D y ) of the VBS order 
parameter. The columnar VBS operator for x-direction 
bonds are defined as 



(7) 



and D y is defined analogously. An SSE-sampled con- 
figuration can be assigned definite "measured" values 
(D x ,D y ) by the operator-counting procedure discussed 
above in the context of the susceptibility (J3j). We ac- 
cumulate the probability distribution P(D x ,D y ), which 
reflects the nature of the VBS fluctuations. In analogy 
with XY models with dangerously-irrelevant Z4 pertur- 
bations [!o|, one would expect the four- fold symmetric 
VBS distribution to develop signatures of U(l) symme- 
try. This has previously been observed when approaching 
the quantum-critical point at T = 0. We now approach 
this point by following the T > critical line. Fig. [5] 
shows results for several combinations of the system size 
and the coupling ratio. While clearly four- fold symmet- 
ric distributions apply for large the histograms become 
more circular-symmetric as the quantum-critical point is 
approached. As at T = [HI, one would expect the 
distribution to be effectively U(l) symmetric when L (or 
some other the course-graining scale) is less than a lengt- 
scale A, with A 00 as q — >• q c . For the system sizes 
studied, we are below A at q = 5/3, while for the larger 
q in Fig. [5] the system sizes exceed A. These observa- 
tions provides direct evidence for U(l)-symmetric VBS 
fluctuations leading to the large v found here close to q c . 

Discussion — All our calculations show consistently 
that the thermal VBS transition in the J-Q3 model has 
critical exponents varying in a range expected in a par- 
ticular subclass of c = 1 CFTs. The exponent n is con- 
stant at r] = 1/4, in agreement with weak universality, 
and v grows rapidly as the quantum-critical point is ap- 
proached, indicating an emergent U(l) symmetry of the 
VBS order parameter and a KT transition obtaining in 
the limit T c — >> + . While we cannot strictly rule out 
a change of behavior to a first-order transition for very 
low temperatures 11, 14, [22|, there are no indications 
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FIG. 5. (Color online) Dimer-order distribution P(D X , D y ) 
for system size L = 32 (left panels) and L = 64 (right panels) 
in the close vicinity of T c . The coupling ratios (temperatures) 
are q = 10 (T = 0.29) in (a),(b); q = 10/3 (T = 0.218) (c),(d); 
q = 5/3 (T = 0.08) in (e),(f). In (f) the distributions is ef- 
fected by unequal sampling (due to long QMC autocorrelation 
times) in different angular sectors. 

of this in any of our results. Note, in particular, that in 
finite-size scaling at a first-order transition one should see 
v = 1/d [41] . where d is the dimensionality (i.e., d = 2 
in our case when T c > 0). Instead, at the lowest T c 
reached here, v w 3. We expect that the same behavior 
should apply also in the J-Q2 model, but that cross-over 
behaviors associated with the proximity to the quantum- 
critical point for all Q2/J in that model may make it 
difficult to extract the exponents there [34[. 

The significance of establishing the nature of the T > 
critical line is that it puts the phase diagram of the J- 
Q model firmly within an established CFT. In the limit 
T — >> + , the effective (2 + l)-dimensional system, ob- 
tained in a quantum-classical mapping through the path 
integral, can still be considered finite in the "time" di- 
mensions, and, thus, the KT scenario can apply. Exactly 
at T = the effective system is fully 3D and a different 
criticality must apply (that of the putative DQC point). 
Since microscopic details should not matter, by univer- 
sality our results should apply to a wide range of VBSs. 

The non-commutability of the limits L — >• 00 and 
1/T — >• 00 is also associated with interesting cross-overs, 
which we have observed here but not studied in detail. 
Further investigations of this aspect of the AFM-VBS 
transition are warranted. 
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